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Abstract 

We study the dynamic and thermodynamic origin of non-canonical equilib- 
ria, and we discuss their connection with the generalized central limit theorem 
and the micro-canonical Boltzmann principle. We reach the conclusion that 
the zeroth law of thermodynamics and the Boltzmann principle are fulfilled 
thanks to an apparent fault turned into a benefit: the dynamic approach can 

only produce a truncated form of inverse power law equilibrium. 
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I. INTRODUCTION 



This paper aims at establishing the most convenient theoretical framework accounting 
for the emergence of non-canonical form of equilibrium. As pointed out by the authors of 
Ref. one of the big merits of the non-extensive thermodynamics of Tsallis is that of 
making popular the discussion of non-canonical equilibria that a few years earlier might have 
caused the severe criticism of the referees We address the same problem from within 
the dynamic perspective of Ref. |^ and we reach conclusions that support theoretically the 
existence a form of non-canonical equilibrium, in apparent accordance with the reasons of 
the success of non-extensive thermodynamics [Q. Actually, as we shall see, our conclusions 
significantly depart from the tenets of non-extensive thermodynamics d la Tsallis, and 
restate the importance of the ordinary Boltzmann principle 1^-0], in accordance with the 
point of view of Gross, illustrated in one of the papers of these Proceedings [Q. 



II. JAYNES APPROACHES TO LEVY PROCESESS 

It has been recently pointed out that the adoption of the method of entropy maxi- 
mization d la Jaynes |T(],TT|, with the Shannon entropy replaced by the Tsallis entropy, does 



not yield directly the Levy distribution, but a probability density function n(a;) that for 
reiterated application of the convolution yields the stable Levy distribution. Here we show 
that, in principle, the method of entropy maximization yields exactly the Levy distribution, 
if a deliberate use is made of the informationa nature of the Jaynes method ||lO| , p!T| . 

Let us discuss this way of proceeding in the Gaussian case, first. Let us imagine that the 
problem to solve has to do with determining the most probable form of a square summable 
function f{x) belonging to the real axis [—00,00]. We know that this function is symmetric 
around x = and we know the first two non-vanishing terms of its Taylor series expansion 
about X = 0, 

f{x) = ^-^x^ + .... (1) 



It is convenient to stress that in this case the information available to us is expressed by 

/(O) = ^ (2) 

and 

It is apparently difficult to proceed with the method of entropy maximization to guess 
the unknown form of this function. However, this is made easy, if we move from the x- 
representation to the k- or Fourier-representation of the function /. We notice, in other 
words, that the information available to us can be expressed under the form of moment 
constraints if we quit the representation of / as a function of x, and we focus on the Fourier 
transform of /(x), denoted by f{k). In fact, the knowledge of cq and C2 can be used to set 
the constraints 



dkf{k)dk = Co (4) 

and 

' dkf{k)k^dk = C2. (5) 



According to the principle of entropy maximization [|^, we have to look for the maximum 
of the Shannon entropy 

S[f] = - j dkf{k)lnf{k), (6) 

while taking into account the constraints of Eq. and Eq.(|^) by means of the Lagrange 
multiplier method. This yields 

f{k) = Ae-^'\ (7) 
this resulting form being the Fourier transform of an ordinary Gaussian distribution. 
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The derivation of the Levy distribution from the proper extension of these arguments 
is easy. We assume that the second piece of information available to us, rather than being 
expressed in terms of the second-order derivative of Eq. is given by 

r/(x)U=o = -^, (8) 



while the form of the first piece of information of Eq. is kept unchanged. Note that the 
symbol -^^^ denotes the symmetric fractional derivative |jl3|, defined by its action on the 
Fourier space, which reads: 

-^|^/(x)} = -ifcr/(fc). (9) 

This means that in the Fourier space the constraint on the second moment of Eq. is now 
replaced by 

J dkf{k)\k\''dk = Ca. (10) 

In this case the method of entropy maximization yields 

/>) oc e--l^^ (11) 



which is well known to be the Fourier transform of an a-stable Levy process ||14|| . 

In conclusion, we have solved the problem raised by Shlesinger and MontroU [0] about 
the derivation of Levy processes from a maximum entropy formalism, without departing from 
the extensive form of Shannon entropy and without using strange logarithmic constraints, 
as proposed by these authors However, this interesting conclusion does not establish 



the thermodynamic character of Levy processes, since the information approach in this case 
might depart from the Boltzmann principle 



III. ON THE WEST-SESHADRI STATISTICS 



To set the dynamic approach in the proper perspective, we remind the reader about the 
deep connection between the Boltzmann principle and dynamics established by the authors 
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of Ref. 0. An oscillator, with coordinate x and velocity v is coupled to a dynamic system, 
called booster to emphasize the fact that no thermodynamic property is already used, as 
it happens with the ordinary thermal baths. The Hamiltonian coupling is given by nxC,, ^ 
being the coordinate of one of the particles of the booster, referred to as doorway variable. 
The booster is a dynamic system in a condition of strong chaos. We do not assign to it 
any thermodynamic property, but a given energy E and we assume that the condition of 
strong chaos makes it reach the micro-canonical state. We write the Liouville equation of 
the whole system, oscillator plus bath, and we derive from this equation, with a projection 
method, the equation of motion of the oscillator. We prove that, under the condition of time 
scale separation, this reduced Liouville equation becomes equivalent to the Fokker-Planck 
equation, thereby leading to a canonical form of equilibrium. No use of thermodynamics 
has been done to reach this important conclusion, but only dynamic properties have been 
invoked 0. At this stage we assume that the width of the resulting oscillator equilibrium 
distribution can be identified with the temperature T of the booster, and we obtain 

ksT = [—lnW{E)—ln <e> +Re^i{u)]-\ (12) 

where W{E) denotes the volume of booster multidimensional surface with energy E and 
$^(co') is the Laplace transform of the correlation function of the doorway variable ^, eval- 
uated at uj, the oscillator frequency. We know that W{E) is proportional to E^ , N being 
the number of degrees of freedom of the booster. Thus, in the thermodynamic limit only 
the first term within the square bracket of Eq. (|12]) survives, thereby recovering the ordinary 
form of Boltzmann principle. It is evident that this nice connection between dynamics and 
thermodynamics is made possible by the fact that, using the language of the advocates of 
non-extensive statistical mechanics 0,^, extensive conditions apply. This is so because of 
the time scale separation between system of interest and booster, resting on the fact that the 
function ^^(t) is integrable. Furthermore, the oscillator of interest is coupled to only one par- 
ticle of the system. The coupling with all the particles of the booster would prevent us from 
recovering the ordinary form of Boltzmann principle. In this paper, we focus our attention 
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on the case where the former extensive condition is broken, and hm(_,. ^^{t) — cons t/t^, with 
P < 1. We shall see that this dynamic condition generates a special kind of non-canonical 
equilibrium. To prepare the ground for this interesting form of non-canonical equilibrium, 
illustrated in Section II B, we have to discuss first the case of free diffusion. 

A. Free Diffusion 

We have now to address a crucial issue, that concerning the subtle difference between 
Levy flight and Levy walk. Let us imagine a random walker that walks according to the fol- 
lowing prescription. At regular intervals of time, 0, T, 2T, . . ., we draw the random numbers 
^1, ^2, ^3, • • •, playing the role of random velocities. These random numbers are characterized 
by the probability density P{r)) given by 

1 

the factor of 1/2 taking into account the fact that the probabihty for the walker to make 
jumps in the positive direction is the same as that of making jumps in the negative direction. 
Note that our discussion rests on the special case where the first moment of this distribution 
is finite, while the second is infinite, namely, 

2<n<3. (14) 

The reason for this choice is transparent. The fact that < 3 ensures that the second 
moment of the distribution is infinite, thereby making it possible for us to depart from the 
Gaussian basin of attraction, which would correspond to the case of canonical distribution. 
On the other hand, we have to ensure that the first moment is finite, so as to establish, as 
we shall see, a connection with a satisfactory dynamic approach. 

Let us consider the case where the random drawing is carried out a given number of 
times, n. This means that the random walker, moving along an one-dimensional path, the 
X-axis, at "time" n is found in the position 
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x{n) = {r]i+ri2 + . . . r]n)T. (15) 

Let us imagine now that this process is repeated an infinitely large number of times, so as 
to build up a distribution density yielding the probability for us to find the random walker 
at the position x at time n, p{x, n). According to the generalized central limit theorem ||15|| , 



the Fourier transform, p{k, n), in the limiting condition of very large n's, gets the analytical 
form 

p{k,n) = exp{—b\k\"n), (16) 

where 

a = /i — 1. (17) 
In the physical condition corresponding to Eg . ([T^) , the diffusion coefficient b reads 

It is important to point out that the asymptotic form of Eq. (|16D is reached after applying 
the convolution procedure for a finite number of times, ricrit, which is fixed to be of the 
order of 10 [|l^. Thus, we make the assumption of considering a coarse grained time scale 
n, where also the infinitesimal change d-h fits the important property dn » Ucrit- 

Let us now consider another random walk prescription. This has to do with drawing the 
random numbers r^'s, with the probability density ipi^r) given by 

^(^) = (/^-l)(^^- (19) 

This means that we can build up an infinite sequence of numbers, which are then used to 
make a random walker walk with the following rules. At time t = 0, when the first random 
number, ti, is selected, we also toss a coin to decide whether the random walker has to move 
in the positive or in the negative direction. The random walker walks with a velocity of 
constant intensity W. Thus, tossing a coin serves the purpose of establishing whether the 
velocity of the random walker is W (head) or —W (tail). This condition of uniform motion 



lasts for a time interval of duration ri. At the end of this condition of uniform motion, a 
new number, T2, is randomly drawn, and a new velocity direction is established by another 
coin tossing. It is important to stress that the physical condition of Eq.(|T^) corresponds to 
the non- vanishing mean time < r >, whose explicit expression is: 

We denote as event the joint process of random drawing of a number and of coin tossing. 
We consider a time scale characterized by the property t >>< r >. It is evident that the 
number of events that occurred prior to a given time t is given by 

n = (21) 

< r > 

Note that we set also the condition that n > ricrit, thereby implying that rather than the 
absolute time t we are considering the coarse time i defined by 

t = n<T > . (22) 

If we adopt this coarse-grained time scale the position occupied by the the particle 
according to Eq. (|T3p becomes indistinguishable from the random walker position prescribed 

by 

=^in+^2r2 + ...^uTn, (23) 

where denotes a stochastic variable with only two passible values, either W or —W, a 
variable that keeps the same sign for the whole time duration of a time interval between two 
nearest-neighbor events. We can also connect P{ri) to ipi^r) as follows 

rji I I rri 

P(„) = ^(^). (24) 

As pointed out in Ref. 0, this has the effect of making Levy diffusion compatible with a 
dynamic and Hamiltonian derivation. 
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B. Fluctuation-dissipation without a finite time scale 



In 1982 West and Seshadri made an interesting proposal to derive a form of non- 
canonical equilibrium. This is described by the Langevin equation 

^ = -Mt)+VL{t). (25) 



Note that West and Seshadri assumed the variable 77 to be a Levy stochastic process. 
This means that the continuous time t of their treatment must be identified with the coarse- 
grained time of Section IIA. This is an important aspect that has fundamental consequences 
on the dynamic realization of Levy processes, and, consequently, on our dynamic approach 
to non-canonical equilibrium. The equilibrium distribution emerging from Eq. ( P^D can be 
easily evaluated by noticing that the probability distribution p{x, t) is driven by the following 
equation of motion 

This is the density picture equivalent to Eq. (^). Both equations afford an attractive 
picture of dynamics driven by both dissipation and fluctuation. In the ordinary case the 
fluctuation process is described by a second-order differential operator. The anomalous 
case here under discussion forces us to replace the second-order derivative with a fractional 
differential operator. The Fourier transform of Eq. (p6[) obeys the time evolution equation 

^^p{k, t) = -b\krp{k, t) - ik-^PiK t). (27) 

This equation yields the equilibrium distribution 

f){k,oo) = exp\ — ^|A;|"), (28) 
V "7 / 

which we refer to as West-Seshadri (WS) non-canonical equilibrium. It is important to point 
out that this form of equilibrium distribution has slow tails inversely proportional to \x\^. 
In other words, in the case 2 < < 3 the equilibrium distribution keeps unchanged the 
power law nature of the original fluctuation. 



At this stage we have to fit the request of making our treatment compatible with a 
Hamiltonian derivation . As pointed out in Ref . , in accordance with the Hamiltonian 



formulation advocated by Zaslavsky |[T^, this condition is fulfilled by replacing Eg. (p5D with 

^ = -7x(t)+e(t), (29) 

where ^(t) is the stochastic process described in Section IIA. This has apparently the effect 
of making the WS statistics compatible with a Hamiltonian derivation. However, Eq.(p9|) 
yields an equilibrium that is not exactly equivalent to the predictions of the WS statistics. 
In fact, it is straightforward to prove that the trajectories moving always in the same direc- 
tion, namely, those exploring the largest distances from the origin, cannot overcome, in the 
positive direction, the distance x = Xmax = — and, in the negative direction, the distance 



w 



. In other words, in the long-time limit, - >> T, we realize a truncated 



Levy equilibrium. In Section V we shall see that this property becomes the key ingredient 
to ensure thermalization between a Gauss and a Levy system. 

IV. THE MICRO-CANONICAL BOLTZMANN PRINCIPLE AND THE 
GENERALIZED CENTRAL LIMIT THEOREM 

The central idea of this section is borrowed from Rajagopal and Abe [|l9l. In fact, these 



authors made the remarkable observation that according to Khinchin ordinary statistical 



mechanics rests on the central limit theorem thereby making it reasonable to expect that 
the non-extensive statistical mechanics is based on the Generalized Central Limit Theorem 
(GCLT) |15|. We find this observation very appealing and we want to discuss in this section 
its consequences from within our perspective that, as mentioned in Section I, is also based 
on the micro-canonical Boltzmann principle, in accordance with other authors 

Let us assume that we know the energy distribution of a small subsystem of a macro- 
scopic system that is assumed to obey the micro-canonical Boltzmann principle. The small 
subsystem is denoted, for the sake of simplicity as particle. Let us denote by p{e)de the 
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probability that the particle energy e is found in the small interval [e,e + de\. We do not 
take position on the explicit form of this energy distribution. We only make the assumption 
that 

7- / \ const 

lime^ooP{e) = (30) 

with < 1/ < 2. This is compatible with both Leyy and TsaUis statistics. Let us assume 
that the macroscopic system consists of N independent particles and let us define the energy 
per particle, e: 



(31) 



N 

Since the N particles are independent the ones from the others, we obtain for e the following 
expression 

PN{e) = / m^3P{e,)]5 [e - Si j . (32) 

The characteristic function of the distribution Pjv(e)), PN{k), is related to the characteristic 
functions of the single particle probability distribution p{ej), p(efe), by 

PN{k))^p''{^). (33) 

The key aspect of the search we are doing is the following one. We study the limiting case of 
N ^ 00 ior the purpose of assessing if the characteristic function becomes equivalent to the 
Fourier transform of a delta of Dirac. If this happens, then the non-canonical equilibrium 
under study is compatible with the micro-canonical principle. If it does not, the non- 
canonical equilibrium distribution is found to be incompatible with the micro-canonical 
principle. 

We skip the discussion of the case 1/ < 1, which is proved to be incompatible with the 
micro-canonical condition, and we focus our attention on the condition: 

l<u<2. (34) 

In this case the first moment is finite and is denoted by the symbol a, namely 
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a =< e >= 



j ep{e)de 



(35) 



and the Levy-Gnedenko theorem [jl5[ affords an analytical expression for the asymptotic 
distribution of the variable 

Ef=i efc - Na 



u 



N- 



(36) 



We denote by Ljq{k) the characteristic function of the variable u, and we obtain 

limjq^aoLNik) = exp 



ik-i-h\kY { 1 + /5-5^tan(— ) 
\ \k\ 2 



(37) 



and 



LN{k)) = e" 



1 1 k 



(38) 



thereby yielding 



limj^^ooP^ i — r) = exp ik'y + ikaN^^'^ — b\k\'^ { 1 + i(3— tan — 
N- [ \ \k\ 2 ) _ 



and consequently: 



limN^^p'' [—) = exp 



ik-fN--^ + ika - h\k\''N^~'' I 1 + ip— tan 



vrz/ 



(39) 



In this case we see that increasing rather than a broader and broader distribution makes 
the righ hand side of Eq.(^) identical to exp{ika), namely the Fourier transform of a delta 
Dirac. Thus, this case is compatible with the micro-canonical equilibrium. 

At this stage it would be straigthforward to prove that the Tsallis non-canonical equilib- 
rium is compatible with the Boltzmann principle provided that the entropic index q fulfills 
the condition 1 < g < 7/5. The WS statistics, on the contrary, would be incompatible with 
the Boltzmann principle. However, as we shall see in Section V, the dynamic approach to 
WS statistics, thanks to the fact that the long tails are truncated, fits both the Boltzmann 
principle and the zeroth law of thermodjTiamics at the same time. 
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V. THERMAL EQUILIBRIUM BETWEEN A NON-CANONICAL AND A 
CANONICAL SYSTEM AND CONCLUSIONS 



All the moments of a truncated Levy process are finite, thereby fitting the Khinchin 
prescriptions for an equilibrium distribution to be compatible with the Boltzmann principle. 
The same property makes this form of non-canonical equilibrium compatible with the zeroth 
principle of thermodynamics. In Ref. [|I|] it has been shown that the finite second moment 
of the non-canonical distribution can be evaluated analytically and its explicit expression is 

< x'ioo) >=< e > T^-'exp[^{^^ - m ^^^~,%^^\ (40) 

where T{a, z) is the incomplete Gamma function. This conclusion can be turned into a 
benefit. We think in fact that it makes it possible to establish in a natural way the thermal 
equilibrium between a system with Levy statistics and one with ordinary Gauss statistics. 



This naturally emerges from the theoretical approach established years ago in Ref. ||2T| to 
study the process of heat tranfer from a hotter to a warmer system. Let us imagine that 
a Gaussian oscillator, with temperature T, is weakly coupled to a truncated Levy process, 
whose temperature has to be assessed using the Gauss system as a thermometer. The 



thermometer equilibrium is shown ||2T| to depend only on the second moment of the Levy 
system, thereby making it possible to establish naturally the thermalization between the two 
systems, and also to measure the Levy temperature by means of the Gauss thermometer. 
In conclusion, the adoption of the dynamic perspective of Ref. , extended to the case of a 
booster producing fluctuations with infinite correlation time, yields a non-canonical form of 
equilibrium, which is compatible with the Boltzmann principle and with the zeroth principle 
of thermodynamics. 
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